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Abstract On the basis of the study of inclusive hadronic events, two methods are 
adopted to determine the number of produced 7(25') events collected by BES in 2001—2002 
run, which is 14.0 x 10° with the uncertainty of 4%. 
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1 Introduction 


The BEijing Spectrometer (BES) is a general purpose solenoidal detectorl!] running at Beijing 
Electron Positron Collider (BEPC). The beam energy of BEPC is in the range from 1.5 GeV to 
2.8 GeV with a design luminosity of 1.7 x 10°'cm7~?s7! at 5.6 GeV center of mass energy. The 
main physics goal is to study the charm and 7 physics. During 2001-2002 years’ running, about 
14 million 7(2S) online hadronic events have been collected!. On the basis of this large data 
sample, many physics analyses could be performed with an unprecedented precision. 

The determination of the offline total number of ~(2S) event, N, JOS) is a foundational work 
in physics analysis, and in turn is the foundation of the further analysis study. In Y(25) physics 
analysis, the calculation of the absolute branching ratio depends on N, pS)’ whose error will be 
directly accounted into the error of the branching ratio of any being studied channel. Therefore, 


it is essential to work out the Noes) accurately and reliably. In principle, any decay channel 
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'Totally 3000 runs (RUN20050-RUN23085) are taken, among which only those with Iquality = 2,3 are 
used to determine the total number of the (2.5) events. Here Iguatizy denotes the run quality and value 
2 and 3 indicate that the run quality is fairly reliable. 


with known branching ratio could be used to evaluate the total number of w(2S): 


b 
TOT _ A 
Y(25) ep- Bp’ 


where Nps is observed number of final state f for a certain decay channel, ep and By are the 
corresponding efficiency and the branching ratio. It is obvious that the larger the branching 
ratio and the smaller the corresponding error, the more reliable the total number is. On such an 
extent, the inclusive hadron final state is a favorable process for the total number determination. 
The only disadvantage here lies in the difficulty to eliminate all kinds of backgrounds throughly. 
Therefore the meticulous studies have been made for the hadron event selection. 

In the following sections, the hadron event selection is discussed firstly, then two methods 
are utilized to determine the total number of 7(2S') and the uncertainties from various sources 


are studied. At last, the final result is given. 


For clearness and convenience, some notations which are to be used afterwards, are listed in 
the Table 1. 


Table 1: Notations 


Symbol Meaning Superscript | Meaning Subscript Meaning 
N Observed Number T Total h hadron final state 
m Selected Number P Peak region e ete” final state 
n “Pure” Hadron Number R Resonance 
€ Efficiency C Continuum 
o Cross section 
L Integrated Luminosity 


In addition, there are two elementary relations among the five quantities N, n, €, o, and L, 
that is 


or N=L-o , (1) 


or n=WN-e. (2) 


The symbol with a tilde on it (e.g. ñ, N, etc.), denotes the events obtained at the continuum 
region (Eyecam = 3.665 GeV), while others denote the events obtained at the resonance region 
(Eveam = 3.686 GeV). 


There are also two frequently used equalities, the first one for variables of the same process 
at different energy points: 


N/a L n/(é-G) 
N/o L nf(e-o) ’ 


s/s 
œ: 
Q: 


= —— (forč=e) ; (3) 


Event 


the second one for variables of different process at the same energy points: 


wie N?! Lo -of w gagi 


nifel NI Lool oF’ Rt Told | (4) 


2 Hadron Event Selection 


For the hadron event selection, the detail information could be found in Refs. [2] and [3]. There 
is no particular event topology to require; instead cuts are made to reject major backgrounds: 
cosmic rays, beam associated backgrounds, two-photon process (y*y*), mis-identified “hadron” 
event from QED processes of ete~ — IT1~,1 = e, u, T, and ete” — yy followed by y conversion, 
and so forth. Most of these kinds of event have salient topology and could be eliminated by 
proper criteria. Events with at least two well reconstructed charged tracks within | cos 8| < 0.8 
are selected (that is Ngooa > 2). The total energy deposited by an event in the BSC (Esum) 
is required to be larger than 0.36 Ebeam, in order to suppress the contamination from two- 
photon processes and beam associated backgrounds. Events with all tracks pointing to the same 
hemisphere in at least one of axial directions (x or y or z direction) are removed to suppress 
beam associated backgrounds. (This requirement could be expressed quantitatively as Issi > 1, 
where Issi is called the squared spatial distribution index.) For two-prong events, two additional 
cuts are applied to eliminate possible lepton pair backgrounds. The number of photons must be 
greater than one (that is Neal y > 2), and the acollinearity between two charged tracks, a Acol, 
must be greater than 10 degrees. 

After the event selection, the fitting of double Gaussian plus a polynomial is applied to 


eliminate the remained background from beam associated backgrounds, see Fig. 1 and 2. 
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Figure 1:  Vertex-Fit distribution Figure 2: Vertex-Fit distribution 
(Eem = 3.686 GeV). (Eem = 3.665 GeV). 


The global error analysis approach is adopted to obtain the uncertainties of selection cuts [3], 
as listed in Table 2: 


Table 2: Error of hadron event selection. 


Requirement Error 
Ngood 2 2 2.765 % 
Esum > 0.36 Ebeam 2.544% 
Issi <1 0.173 % 


Nrealy > 2 (for Ngood =9 ) 0.042 % 
Cte Z 10° (for Ngoa = 2 ) 0.157 % 
Sum 3.765 % 


In fact, there are many processes that could lead to hadron final state at ~(2S) resonance 


region, they could be divided into seven categories: 


ete” — hadron (CH) , (5) 
ete” — (2S) — hadron (RH) , (6) 
ete” — J/w—hadron (J/-H) , (7) 

(8) 

(9) 
) 
) 


ete” — Ttr — hadron (7-CH) , 


ete” — (2S) —7tr7 — hadron (7-RH) , 9 
ete” — hadron* (CH*) , (10 
ete” — w(2S) — hadron* (RH*) , (11 


where C represents the continuum process, R the resonance process, H hadron event, and H* 
indicates the event which survives all aforementioned hadron selection cuts and is left in hadron 
sample from processes (10) and (11). Among above seven categories, only the hadron event 
from the first two processes are “pure” hadron event at 7(2) peak region while others should 
be treated as backgrounds. Since hadron event from different process has almost the same event 
topology, the theoretical estimation method is used to evaluate the contamination of such kinds 
of backgrounds. According to the analysis in Ref. [3], two factors are introduced to subtract 
the hadronic background. If the pure hadron number is denoted as n and the selected hadron 


number denoted as m, then it could be obtained 


YOM rad = Nad , (12) 
RM a = neg , (13) 
where 
yo = l- fay — fo = 0.855 , 
yr = 1-fr=0.999 , 
with es aie 
fry = Thad ` Ehad 
oad ` fad 
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C.C 
fo = a 
= C C 
k=r,e,u,y,y* Chad ` Chad 
25 
g g 
fr = R 


p(25) R 
l=r,e,uy Bh Eh 


Here all efficiencies are obtained from Monte Carlo simulation [3], and symbols 7, e, u, y and * 
denote final states T*T~,eTe~, u" u, yy and Ņ*9*, respectively. For the continuum process, 


the production cross section of could be obtained from the corresponding Monte Carlo gener- 


ator. For Fe and ollt, they could be calculated by theoretical formulae with corresponding 
resonance parameters obtained from the scan experiment data [2, 4], which is about 15 nb and 
1 nb, respectively. Since Monte Carlo generator does not give the cross section for resonance 
process, the ratio of branching fractions is used in the calculation of factor fr [3]. It should 
be pointed out that because the variation of ofl is fairly smooth at (2S) region, refer to 
Fig. 3, the contamination from J/7 decay is suitable to be treated as the continuum-like hadron 


background. 
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Figure 3: Cross section curve of J/z resonance. 


3 Determination of Total Number 


3.1 Principle 


The branching ratio of hadron final state, denoted as Bp, can be acquired from PDG list [5] 
or from BES scan results [2]. If the number of selected hadron event from 7(2S) resonance is 
Ny(28)—had., then 

Myp(2S)—had. 


TOT _ 
Ny23) = e£ : Bp, (14) 
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However, the number of selected hadron event at a certain energy point is the combination of 


two parts (refer to Fig. 4), one from resonance process and the other from continuum one, that 


is? 


Nhad = My(28)shad. F Nete--shad. 5 
4 4 4 (15) 
rm = nR) + nk 


Therefore the key issue here is to distinguish the nit from the ne. There are two methods, the 
fraction subtraction method and the normalization subtraction method, can be used to figure 


out the number of resonance event. 


Figure 4: Fractional sketch of hadron number of different process. (The cross-hatched 


area in the figure indicates the data taking region.) 


3.2 Fraction subtraction 


Refer to Fig. 4, if the ratio of ng to n could be estimated: 


C 

f= 

=o 
h 


together with Eq. (15), the n? can be calculated as 


P ng 
Nh 1+f (16) 


As to the factor f, from Eq. (4), it is easy to acquire 


nl €. oC E. oC 

h _ &h PR ah ye ae = 
nE a h UR gk "h >? 
h h` Oh h` Oh 


2 Because of energy spread and shift, the data are actually taken within a region rather than at a 
fixed point (notice the cross-hatched area in Fig. 4). So “P” instead of “R” is adopted in the following 
formulae in order to denote such an effect. 


so f can be expressed as 


f= R -R * (17) 


Combing Eqs. (14) and (16), it can be obtained 


T 
NIOD = = . 18 
2S) B,- ef. (14 f) (18) 


3.3 Normalization subtraction 


The data taken far from the resonance region could be treated as the data of continuum?, refer 


to Fig. 4, from Eq. (3), the number at continuum region ñF could be transformed into that at 


resonance ng by a luminosity normalization factor 


c -C 
nk = TÙR , fr= 


where fr is the transformed factor. Combining with the relation n? = n? + n@, the resonance 
h h h 


number could be worked out 
P T C T A 
Ny, =, — ÍT AR = Nh, ‘(1 — fr- oF) 

h 


In the expression of fr, the luminosity L is usually calculated by the continuum ete~ event (i.e. 
Bhabha event): 


Similar to Eq. (16), at the continuum region, the number of events from the continuum process 


is expressed by the number of total selected events: 


T oat >T 


tag? PG n 
with 
p=% ii (19) 
j= 5-2 , (20) 
hei. A 


It should be noticed that the relation e = ERC has been usedź in the above calculation. 


3Strictly speaking, within the scan range, any data are from two processes, resonance and continuum. 
The effect of resonance to continuum could be taken into account by factor fe, fe, fn in Eqs. (19), (20), 


and (21). 
4The relation a a a is exact for the ete” final state, whose event selection cuts are energy 


independent; but for the hadron final state, the relation is only an approximation. 


T 


Put all together, the total number can be calculated as 


T 
TOT _ Ph 

Nps) = B, -2 (1— Fr: Me) , (22) 

where 7 2 
Fp = One U+), 1 (23) 

GE-of (1+ fe) (14+ fn)’ 
T >T 

n n 
Mp = eh . 24 
eh nT nT ( ) 


with factors fe, fe, and fa are given in Eqs. (19), (20) and (21). 


3.4 Correction 


As mentioned in section 2, after the hadron event selection, the selected hadron number m 
instead of pure one n is obtained. The relation between m and n is given in Eq. (12) and (13), 
that is 


yom, me) = (nx , nE) ; 
~ R >R 
yr(mÈ, inj’) = (np > Mp, 
Notice that 
C 
Mp =S Mp EME , (mp, =™Mp) , 
then P Cc P P 
n n n n 1 
mE = 2h 4 Mh PEL (44) , 
YR VC: YR YC YR YC 
that is m 
no h 
1 
pues 
YR YC 
So the formula for the fraction subtraction method now becomes 
T 
NTOT — Mh (25) 


(2S) Br- R- (1+ f) óy i 
where 6,, is defined as 
1 f 
m 
i, = O, 26 
For the normalization subtraction method, the corresponding corrected formula could be 


obtained similarly and the final result is 


NTOT — ip (1 Fr:Men) , 27 

v0) = Tog, ek LIT Men) (27) 
where F 
n m 

Me = =F > (28) 
nt m? 


and 
Fr = Fr/õy ’ y = — ’ (29) 


here Fy is just as that defined in Eq. (23). 


3.5 Numerical Calculation 


By the virtue of Eqs. (25) and (27), N JOS) could be worked out. For convenience, all numbers 


relevant to total number calculation are summarized in Table 3. 


Table 3: Event number, efficiency and cross section. 


Final state hadron ere 
Region Resonance Continuum Resonance Continuum 
Event number mZ = 10634586.0 mz = 13937.7 |n? =1190613 nt = 54415 
mi (Amf”) (Am? = 11688.2) (Am2 = 817.2) 
Efficiency | Resonance eR = 0.753 ER wy eh e£ = 0.761 cel 
process (Aef = 0.012) 
Continuum é = 0.716 EO ee ef = 0.708 erce 
process (AeÈ = 0.010) 
Cross Resonance oÈ = 676.277 Gj, = 0.266 of = 4.053 a = 0.00158 
section process 
(nb) Continuum of = 15.495 o¢ = 15.673 | of =80.428 of = 81.349 
process 
Trigger efficiency { Correction factor 
e9 = 0.99924 Yr = 0.999 Yo = 0.855 


f: the trigger efficiency is obtained from Ref. [6] written by Dr. Fu ChengDong. 


NGOs) is worked out to be either 14.05x10° for Fraction method or 14.04x10 for Normal- 


ization method. 


4 Error Analysis 


4.1 Classification 


Formally, by the virtue of Eqs. (18) and (22), the formula to calculate the total number can be 
written as 


NTOT = Oh, Gi, (30) 
v(25) Bn- eÈ 
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where G; is a correction factor defined as 


1 
G; = —— , for the fraction method; 
G; = 1+ f 


Gə = (1 — Fr- Men) , for the normalization method. 


The error of NGOs) comes from the components of Eq. (30), such as nf, eft. Bp, and G;, the 


error of which will be discussed one by one. 


4.2 Uncertainty of selected number m$ 
For the selected number? mi, there are three sources of uncertainty®: 


1. Fitting uncertainty 
The uncertainty of fitted number could be obtained from the corresponding error of fitting 


parameters which are used to calculate the number (refer to Table 3), that is 


2. Statistic uncertainty 


According to statistic principle, 


1 
Vsta(™Mp, ) = 
Vm, 


3. Selection uncertainty 
According to the study of section 1, the selection uncertainty reflects the inconsistency 


between data (m?) and Monte Carlo (ef), so the uncertainty of e} is also included in this 


T 
V sel (=) = 3.77% 


h 


term which is 


4. Effect due to beam energy fluctuation 
As it is mentioned in footnote 2, the data are actually taken within an enegy range (refer 
to the sketch description in Fig. 4), so the ratio between g / an will vary with the actual 
beam energy which may be different for different beam-injection. Usually, each beam- 
injection includes 3—5 runs. As an estimation, all runs are grouped with every 3, 4 or 5 
runs, then the ratios of selected hadron number (n}) to that of ete~ number (n7) are 
worked out, which denoted as reli) with 7 indicating the grouped run number. Fig. 5 
shows the rpe distribution for different grouped-run number. The maximum value of rp,¢(7) 


corresponds to the peak cross section’. Taking the experimental statistic fluctuation into 


5Hereafter the selected number m instead of pure number n is used in the error analysis, and the 
uncertainty for such substitution is rather small and is to be discussed afterwards. 

°Hereafter the symbol v denotes relative error. 

"Notice 


account, the value 10 is adopted as the position of the peak cross section. The uncertainty 


from the beam energy fluctuation effect is estimated as follows 


vy =< — malt f) 
k mA /L+ FO] 


with 
C C 
F) = -k 


Ep -of Theli) 


For different grouped-run cases, Verp is almost the same, which is 0.23 %. 
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Figure 5: rre distributions for different grouped-run number. 


4.3 Uncertainty of branching ratio Bp 


There are two B, values, one from PDG2002 and the other from w(25') scan experiment: 


B,(PDG2002) 
Br(w(2S) scan) 


(98.10 + 0.30))% , 
(97.79 + 0.15)% 


In the total number calculation the By, value form PDG2002 is used. However as a conser- 
vative estimation, the difference between above two values is used as the uncertainty of Bp, that 
is 


v(Bp) = 0.32% 


then 
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4.4 Uncertainty of correction factor G; 


For G4, the uncertainty due to different f could be calculated as follows 


1 1 
A= Nges (S) — Nges’) E 1+f 1+f _ =F 
Nis) (Pf) ES Lf 


According to the definition of f, Eq. (17), the uncetainty of f mainly comes from the statistics 
of eC and eÈ, both of which are 1/v50000. Notice f = 0.02179 is a small quantity and the 


maximum difference between f and f’ is also small, so 
v(G&i) + |1 — f|- A 0.88% , 
where A = |f — f'|max = 2/v 50000. 


For Go, the uncertainty due to different G2 could be calculated as follows 


Pee ee Nias) (G2) 


Nyas)(Ga) 


Notice Fr(= 0.932) approximates to one, so the difference of total number calculated with 


Fr = 1 and with Fr # 1, is treated as the error from factor Fyr, that is 
Men: (Fr — 

ae | eh: (Fr — 1) 
1— Fr: Men 


Fr - Men — Fh - M! 
v(G2) = -|= a 


1—Fr- Men 


= 0.20% 


Notice Men consists of two pairs of ratio, so the systematic error of numerator and denominator 


will cancel automatically, only the statistic error is left, which is 


l 1 1 l 
AMa= op + 5p + art 7 = 0.95% , 


so the error due to different Men could be calculated as 


Fr - (Men — M!,) | Fr - AMen 


zx 0.91% . 
1— Fr Men a 1 — Fr Men 


(Me) =| 


Then the uncertainty for G2 is 


v(G2) = y v?(Fr) + v? (Men) = 0.94% . 


4.5 Other uncertainty 


The other effects which could lead to the uncertainty include the correction factor yc, yr, trigger 
efficiency and so forth. All these uncertainties are regarded as less than 0.5%. In addition, Monte 
Carlo sample is used as input to test the bias of our method. The error from such bias is about 
0.6%. 


4.6 Summary 


Put all things together, the synthetic uncertainties are summarized in Table 4. 
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Table 4: Error summary. 


Source Uncertainty 
a Fitting 0.11 % 
Statistic 0.03 % 
(&e?*) Selection 3.77 % 
Eveam fluctuation 0.23 % 
Bh, 0.32 % 
G; Gi 0.88 % 
Go 0.94 % 

Other (0.560.6) % 
Total Method 1 3.97 % 
Method 2 3.98 % 


Note: Method 1 for Fraction method; Method 2 for Normalizaion method. 


5 Result and Discussion 


The final total number of w(2S) event with corresponding error is 


ror _ | 14.05 x (143.97%) x 10° , (Fraction Method) ; 
1+3.98%) x 106 , (Normalization Method) . 


Notice Eqs. (18) and (22), two methods, the fraction method and the normalization method, 
correlated closely with each other and the difference between two methods for the central value is 
actually less than one per mille. Furthermore, the difference of the uncertainty for two methods 
is also fairly small. Therefore, the central value of the offline total number of w(2S) event could 


be regarded as 14.0 x 10°, and the uncertainty could conservatively be estimated as 4%. 


Thanks are due to Prof. J. C. Chen and Prof. D. H. Zhang for their help about Monte Carlo 


simulation. 


References 


[1] BAI Jing-Zhi et al. Nucl. Instr. Meth., 1999,A344:319-334; 
BAI Jing-Zhi et al. Nucl. Instr. Meth., 2001,A458:627-637 


[2] BAI Jing-Zhi et al. Phys. Lett., 2002, B550: 24-32 
[3] Mo Xiao-Hu, Study of Inclusive Hadronic Event, (2003.6), BES Memo 


[4] BAI Jing-Zhi et al. Phys. Lett., 1995, B355: 374-380. 


201703.00600v1 


chinaXiv 


[5] Particle Data Group, Hagiwara K et al. , Phys. Rev. 2002, D66:010001 


[6] FU Cheng-Dong, Measurement of the Trigger Efficiency of y’, (Feb. 27, 2003), BES Memo 


14 


